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Abstract 

We derive a generalized Luttinger-Ward expression for the Free energy of a many body system 
involving a constrained Hilbert space. In the large N limit, we are able to explicitly write the 
entropy as a functional of the Green's functions. Using this method we obtain a Luttinger sum 
rule for the Kondo lattice. One of the fascinating aspects of the sum rule, is that it contains two 
components, one describing the heavy electron Fermi surface, the other, a sea of oppositely charged, 
spiuless fermions. In the heavy electron state, this sea of spinless fermions is completely filled and 
the electron Fermi surface expands by one electron per unit cell to compensate the positively 
charged background, forming a "large" Fermi surface. Arbitrarily weak magnetism causes the 
spinless Fermi sea to annihilate with part of the Fermi sea of the conduction electrons, leading to 
a small Fermi surface. Our results thus enable us to show that the Fermi surface volume contracts 
from a large, to a small volume at a quantum critical point. However, the sum rules also permit 
the possible formation of a new phase, sandwiched between the antiferromagnet and the heavy 
electron phase, where the charged spinless fermions develop a true Fermi surface. 

PACS numbers: 72.15.Qm, 73.23.-b, 73.63.Kv, 75.20.Hr 
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I. INTRODUCTION 



Sum rules play a vital role in condensed matter physics. The most famous sum rule - the 
Luttinger sum rule 

HQS 

, defines rigorously the volume of the Fermi surface of a Fermi 
liquid in terms of the density of electrons: 

where vps is the Fermi surface volume, n e the density of electrons per unit cell and D the 
dimension. Historically, sum rules have also played an important role in our understanding 



of strong 
sum rule 



correlated systems. In the context of the Kondo effect for example, the Friedel 

DD 

Ev = ( 2) 

A 

relating the sum of the scattering phase shifts in channels labelled by A to the num- 
ber of bound-states An e helped to establish a rigorous foundation for the Abrikosov-Suhl 
resonance which develops in Anderson and Kondo impurity models. Later, Martin jlo| 

applied the Luttinger sum rule to the Anderson lattice model, to argue that heavy electron 
metals must have a "large Fermi surface" which counts both the conduction electrons and 
also the localized f-electrons. 

1 V/VS ' n e + l. (3) 



(2tt) d 

Today, there is a renewed interest in sum rules, in connection with models of strongly cor- 
related electrons. For instance, in the context of high temperature superconductors, which 
are Mott insulators when undoped, there has been a long-standing debate over whether the 



"large" Fermi surface predicted by Luttinger's sum rule, might be r ep 



aced at low doping, 



by a "small" Fermi surface determined by the number of doped holes [111. Il2j|. Related issues 
arise in the context of heavy electron systems, where experimental advances have made it 
possible to tune through the quantum critical point that separates the heavy electron para- 
magnet from the local moment antiferromagnet. There have been a number of theoretical 
speculations that the Fermi surface of the heavy electron material may jump from "large" 
to "small" at the quantum critical point [ls! Il4 lisj] . Recent ex peri mental work, based on 
de Haas van Alphen measurements |l6[] and Hall measurements[17], provide experimental 
support for this hypothesis, but the idea has lacked rigorous theoretical support. 
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The main difficulty in extending sum rules to strongly correlated systems, is that the 
theoretical machinery used and developed by Luttinger and Ward to derive sum rules applies 
to models with un-projected Hilbert spaces. About eight years ago Affleck and OshikawaQ] 
demonstrated that such sum rules have a more general existence. By using a modification 
of the Leib Mattis theorem, Affleck and Oshikawa showed that the "large Fermi surface" 
which counts both local moments and conduction electrons develops in the one dimensional 
S = 1/2 Kondo model, even though in this case, the ground-state is not a Fermi liquid. More 
recently, Oshikawa jl^ has extended this derivation to higher dimensional Kondo lattices. 
This work suggests that it ought to be possible to extend the Luttinger Ward approach to 
Hamiltonian problems with strong constraints. 

In this paper, we show how the original methods of Luttinger and Ward can indeed be 
extended to strongly correlated models. A key element of this work is the construction of a 
functional relating the Free energy of a strongly interacting system to the Green's functions 
of the projected Hilbert space. Our approach is based on the use of slave particles, such 
as a Schwinger boson description for the local moments in a Kondo lattice. We begin with 
a generalization of the Luttinger Ward Free energy functional, appropriate for systems of 



interacting bosons and fermions 



20, 



which can be compactly written as 



F = TStr [In + EQ] + Y[Q] (4) 

where S_trL4] = Tt[Ab] — TrLAp] denotes the supertrace of a matrix containing both bosonic 
(B) and fermionic (F) components (where the underline notation is used to denote a sum 
over internal frequencies and a trace over the internal quantum numbers of the matrix). 
G = (Go 1 — £) _1 is the matrix describing the fully dressed Green's function of all elementary 
particles and fields entering the Lagrangian, including the slave particles, where £ is the self- 
energy matrix and Go the bare propagator of the fields. The quantity Y[Q] is, diagramatically, 
the sum of all closed- loop two-particle irreducible skeleton Feynman diagrams (Fig. U). 
Variations of Y[Q] with respect to the Green's function G generate the self energy 

SY[G] = -TStr [£53]. (5) 

or 



3 



G b = 
*3x = 



(b) 



(c) 




+ ... 



oc 7; 



FIG. 1: Illustrating the Luttinger Ward functional for the Kondo lattice, (a) Fully renormalized 
propagators for the conduction electrons, x fermions and Schwinger bosons, where a £ [1, N] is the 
spin index and v £ [1, K] the channel index, (b) Interaction vertices showing spin (<r) and channel 
(v) indices and (c) leading skeleton diagrams for the Luttinger Ward functional with dependence on 
jj. The first diagram involves two loops carrying spin and channel quantum numbers, and one pair 
of vertices and is consequently of order 0(KN/N) = O(N). The second diagram has four loops 
carrying internal quantum numbers, and four pairs of vertices, so it is of order 0(N 2 K 2 /N 4 ) = 
0(1). The final diagram has two quantum number loops and three pairs of vertices, and is hence 
of order 0(NK/N 3 ) ~ 0(1/JV). 

where the use of the supertrace as the measure for functional derivatives avoids the need to 
introduce a relative minus sign between Fermi and Bose parts of this expression. Diagram- 
maticaliy, the functional derivative of Y with respect to Q corresponds to "cutting" one of 
its internal lines. If we truncate Y to some order in 1/N, this relationship determines a 
conserving Kadanoff-Baym approximation [2 a] . To leading order in the large iV expansion, 
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the self-energies generated from the functional derivatives of Y are 



5Y 

EbM = -P- 




0(1/N) 



where the cross indicates the line which is eliminated by the functional differential. Each of 
these terms contains a factor 0(1/ N) from the vertices, but the first two self-energies contain 
summations over the internal channel or spin indices, elevating and S x to terms of order 
0(1). In the leading order large N approximation, S c ~ 0(1/ N), so that a consistent large 
N approximation is produced by leaving the conduction electron lines un-dressed. This 
provides an alternative dia gram matic derivation of the Par collet- Georges approach to the 
multi-channel Kondo model 26]. 

Each conserved quantity Q which commutes with the Hamiltonian leads to a gauge in- 
variance of the underlying quantum fields, such that the action and all physical properties 
are invariant under the transformation 

Mt) - e ie( ^c- (8) 

Here, is the gauge charge of the field: it is this quantity which controls the charge of any 
physical excitations associated with the field. Provided this gauge invariance is unbroken 
we show that this leads to a Ward identity at zero temperature 



qG- 



duj 



0. (9) 



where the integral runs along the imaginary axis. Luttinger's original work derived the 
elementary version of this relationship for electrons: here it is seen to hold for interacting 
species of fields. 
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Each of these Ward identities leads to a corresponding, generalized Luttinger sum rule 
given by 



Q = —ImTi 

7T 



(10) 



gin (-0^(0 -z<f)) 

where Gf is the complete Fermionic Green's function, involving all fermions involved in the 
description of the constrained Hilbert space, including those introduced as slave particles. 
In impurity systems, the trace over the logarithms reduces to a sum of phase shifts, and the 
quantity Q is replaced by AQ, the change in Q induced by the impurity In the limit of 
infinite band-width, AQ = 0, which leads to a multi-particle version of Friedel sum rule 

» = (|) (id 

where is the spin degeneracy associated with the trace over internal quantum numbers 
of the field ip^. In lattice systems, the trace over the log becomes the Fermi surface volume 
associated with the field ip^, so that 

vfs(C)" 



C 



q ( n c 



(12) 



(2n) D 

where vfs{C) is the Fermi surface volume associated with field fy. A closely related Fermi 
surface sum rule has recently been derived by Powell et al. for mixtures of fermions and 

n 

bosons in atom traps |27|. 

In this paper we apply these results to the Kondo lattice model. Key to our approach is 
the adoption of Schwinger bosons for the description of quantum spins. Schwinger bosons 
offer a key advantage, because they can describe both the antiferromagnetic and the heavy 
electron ground-states of the Kondo lattice|2j|. The Kondo effect induces a development of 
a retarded interaction in the spin singlet channel between the electrons and the spins. For- 
mally, this interaction manifests itself as the mediating field Xj i n the Hubbard- Stratonovich 



decomposition of the interaction j26t 1291. 13( 

JkSj ■ c] a a a pc jl3 -> (ct a 6 ia ) Xj + -!= (&t a Cj- a ) Xj + ^0 (13) 

Here Sj is the spin at site j and Cj a creates an electron at site j. At long-times the x 
field develops charge dynamics which describe the scaling of the Kondo interaction in both 
frequency and momentum space. In particular, the x~ propagator describes the momentum 
and frequency dependent Kondo interaction, 

-^(k ) w)=fi + E x (k,w)j =^(k,w). (14) 



At each point in momentum space where the Kondo interaction scales to strong coupling, 



in 



ted 



15j to 



Q x develops a pole. In this way, the physics of the Kondo effect is intimately 
the possible emergence of a positively charged, spinless Fermi field - a "holon" 

Under gauge transformations, the conduction and \ fields have opposite gauge charge, 

Qc = 1, q x = -!■ ( 15 ) 

When we apply the Luttinger Ward procedure, we find that the physical charge density is 
given by the total volume of the electron Fermi surface, minus the volume of the holon Fermi 
surface, as follows 

Here, the volume of the electron Fermi surface is given by the conventional Luttinger formula, 
v ?s = - Vlmln[e k + S c (k,0-^)-/i + ^)] = V9(^-£ k ) (17) 

7T z — ' z — ' 

k k 

where /x is the chemical potential, E^ = e k + S c (k, E^) is the renormalized energy of the 
heavy electrons. Now v x is 

x _I^Imln[^ + S x (k,O-^))] = ^0(-J^(k)) (18) 

7T '—^ J K 



V 

77 ' 

k 



where J* K (k) = -<7 x (k, 0-iS). 

The meaning of the renormalized coupling constant in the paramagnetic phase needs 
a little discussion. This quantity describes the residual interaction between the electron 
fluid and any additional spins that are added to the Fermi liquid ground-state. If we add an 
additional Schwinger boson, a given site, the additional spin unit 

remains unscreened, because each channel can only screen one spin unit, and all channels 
are fully screening. If Ub = E(nb + 1) — E(rib) is the energy to adding one additional spin, 
then the residual interaction between the additional Schwinger boson and the conduction 
electrons is J^(ub). Since the additional spin decouples, it follows that Jk{oj) has a zero 
at uj = Ub, J^{uJb) = 0. At higher energies, the residual interaction will become positive, 
ultimately connecting up to the single-ion scaling behavior Jk{oj) ~ 1/\d.(Tk/uj). Since 
passes through zero at u = Ub, it follows that at lower energies uj < u>b, and in particular, 
at zero frequency, the residual coupling must be ferromagnetic, i.e J^(0) < 0. This in turn, 
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means that the Fermi sea of holons is entirely full v x = (2ir) , so that 

v fs 



n. 



1. 



(19) 



(2tt) d 

To preserve the overall charge density, the electron Fermi surface volume is forced to enlarge 
by one unit per spin to "screen" the finite background density of postively charged holons 
(Fig. 121). 

Let us now consider the antiferromagnet. In this case, the Schwinger boson field con- 
denses, and the x fi e ld becomes hybridized with the conduction electrons. In this case, the 
holon Fermi surface annihilates with the expanded part of the heavy electron Fermi surface, 
to reveal a single integrated set of Fermi surface sheets with a "small" volume that encloses 
the total electron count. If the transition from paramagnet to antiferromagnet occurs via 
a single quantum critical point, then our results indicate that the transition in the Fermi 
surface volume at the second order quantum critical point is abrupt. 




FIG. 2: Schematic diagram illustrating the enlargement of the Fermi surface in the Kondo lattice 
(a) coupling of local moments to small Fermi surface (b) each spin-boson "ionizes" into an electron 
and a "holon" , producing a large electron Fermi surface and a filled holon Fermi sea. 



Perhaps the most interesting aspect of this picture, is that it suggests that the background 
field of holons can develop dynamics. One intriguing possibility, is that the holons are 
liberated when the Fermi surface volume contracts at the quantum critical point. Another 
more exotic possibility, allowed by the Ward identities, is the formation of a separate phase 
where the holons develop a Fermi surface, so that v x = x(2ir) D is partially filled. This "spin- 
charge" decoupled phase would be sandwiched between the localized magnetic phase and 
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suggested that the holons might develop 
|. This might provide a mechanism for 



the fully formed paramagnet. Pepin has recent 
a Fermi surface at the quantum critical point . 
such a scenario. A good candidate for this phase is the one-dimensional Haldane S — 1 
chain, coupled to two conducting chains to form a one-dimensional, two-channel Kondo 
lattice. Here, the conjectured intermediate phase would seperate a Haldane gap phase and 
a large Fermi surface phase, leading to Friedel oscillations with wavevector q = 2k^ + 
intermediate between that of the large (x = 1) and small (x = 0) Fermi surface. 

This paper is divided up into various parts. In the first we develop the Luttinger-Ward 
Free energy expansion for the Kondo lattice. Next, we present a new formula for the en- 
tropy of the interacting fluid, expressed exclusively in terms of the Green's functions of the 
particles. We then go on to derive the Luttinger sum rule in a general form for a constrained 
system, which we then apply to the Kondo impurity and lattice model. In the final part, we 
discuss the general implications of our sum rule in the Kondo lattice. 

II. LUTTINGER WARD EXPRESSION FOR THE KONDO LATTICE 

The Luttinger Ward approach was originally based on perturbation theory in the strength 
of the Coulomb interaction. In order to develop these methods for strongly interacting 
systems, we need to take account of the constraints in the theory, requiring a Luttinger 
Ward approach appropriate to gauge theories. In this respect, our efforts to apply the 
Luttinger Ward approach to the Kondo lattice closely parallel recent efforts to extend the 



21] 



the use- 



Luttinger Ward approach to quark-gluon plasmas 

Although the Luttinger Ward functional may be derived non-perturbatively 
fulness of any derived sum rules rests on the direct relationship between the Green's functions 
that enter in the functional, and physical excitations. Strongly interacting systems do not 
generally provide a natural small parameter on which to base a perturbation expansion. 
Our philosophy is that to develop this control, we need to define a family of large-N models 
which, in the infinite N limit provide a mean- field description of the phases of interest. In 
the expansion about this limit, 1/N then provides a small paramater. The large N solutions 
constructed this way, will of course satisfy the sum rules we derive, but the sum rules are 
expected to have a far greater generality. 
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In the context of our interest in the Kondo lattice, we adopt a family of models given by 



Hi 



H = E e k C l ua C kua + ^ E C )uaCjupb) p b ja (20) 

Here, we have adopted the Schwinger boson representation, where bj a creates a Schwinger 
boson at site j, with spin component a G [1, N]. The combination b^ a bj/3 = S a p{j) represents 
the local spin operator and the system is restricted to the physical Hilbert space by requiring 
that 

a 

at each site j. The index v G [1, K] is a channel index that is absent in the physical 577(2) 
model, but which is included to create a family of large N models. The "filled shell" case 
2S = K defines a perfectly screened Kondo lattice, and our discussion will focus on this 
case. Here c- creates a conduction electron of momentum state k and spin a G (1,N). 
c^ a = Yl c iia e ~ lk Xj creates an electron at site j, where M s is the number of sites in the 
lattice. The above model is known to have a controlled large N expansion which contains 
both magnetic and paramagnetic solutions when N is taken to infinity with k = K/N fixed. 

To develop a Luttinger Ward expansion, we shall factorize the interaction in terms of a 
Grassman field Xji 

#/-E-/= [(4Mxl + {b] a c jva)Xjl \ + £ (22) 

jua v ju 

and the Lagrangian becomes 

£ = E c L (9r + ^ c kva + e b) a o T + a,) bja + e 

+ E 7^ [tfiMxl + (b] a c jua ) Xjl ] - E 2SX > ( 23 ) 

Note that the Grassman fields Xj contain no time derivative, and so instantaneously, they 
behave as inert, neutral fields. Here Xj is the static value of the lagrange multiplier used 
to enforce the constraint on each lattice site. At finite temperatures the constraint is only 
enforced on the average, but in the ground-state, the value of the conserved quantities 
will be quantized as a step- function of the Xj. In the range Xj G [A~, A^~], ( where A^ are the 
energy gaps between the physical ground-state and the states with 2S ± 1 Schwinger bosons 
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at site j) the constraint will be precisely satisfied. We shall assume that a translationally 
invariant choice Xj = A is always available. 

For our sum rules, we need to know comparatively little about the nature of the ground- 
state excitations. However, from the the large N solution of this model, we do know that 
there are two main classes of solution to the Kondo lattice 

• Heavy electron phase, where the Schwinger bosons pair-condense in either the Cooper 
channel ((ab ia bj_ a ) 7^ antiferromagnetic interactions) or particle-hole channel 
((cbl a bj a ) 7^ ferromagnetic interactions). In field-theory language, this is a "Higg's 
phase" in which the local U(l) symmetry is broken and both the spin and holon 
excitations can propagate from site to site. 

• Magnetic phase, in which the Schwinger bosons individually condense ((bj) ^ 0) to 
develop long-range magnetic order. In this phase, the Xj fields become hybridized 
with the conduction fields, and Q c and Q x become diagonal members of a single matrix 
propagator. 

In general, we will therefore have to preserve the momentum dependence of the propagators. 
The bare propagators of the theory are diagonal in momentum, and given by 

Gf\k,iv n ) = (iUnTs-Xiy 1 , 

g®(k,iu n ) = -J K (24) 

where we have written the bare boson propagator as a two-dimensional Nambu propagator, 
to allow for the possible introduction of boson pairing, which develops as antiferromagnetic 
correlations begin to grow. Formally, we will combine these three propagators into a single 
propagator given by 

W n T3 - XI 

iu n — e? 



g- 1 



i 

' Jk 



- S (25) 



where E is a matrix self-energy. Following Luttinger and Ward, we will regard g as a 
variational function, and £ as a derived quantity. 

We now consider the effect of tuning up the strength of the Kondo interaction from zero 
to J K , by replacing J K — > aJx where a G [0,1], keeping the chemical potential of the 
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conduction electrons and bosons fixed. Now the partition function is given by 

Z = Tr[e^ H ] = J V[c, b, x ]e~ J » ' Cdr (26) 

where V[c, b, X ] is the measure of the path integral. If we vary a inside this expression, we 
obtain 

so that if F — —T In Z is the free energy, then 

^ = -^B4^> = -(^m], (28) 

where Q x is the Green's function for the X fermion and Tr[(/ X ] = K J2 iuJn % G x (k, iuj n )e luln0+ 
denotes a trace over the frequency and momentum of Q x (here, e tu)n0+ ensures the normal 
ordering of the operators). We next consider the expression 

F = TSt L [\n(-g- 1 )+j:g]+Y[g]. (29) 

We will show how this expression generalizes the Luttinger Ward Free energy functional to 
a mixture of interacting bosons and fermions. Here, we use the notation 

StrL4] = Tr[A B ] - Tr[A F ] = ]T Tr [A B (iu n )}e^ 0+ - ]T Ti{A F {iuo n )}e^ 0+ (30) 

to denote the "supertrace" over the bosonic and fermionic parts of the matrix A. The 
"underline" notation is used to indicate a sum over both the frequency variable of A and a 
trace over the internal quantum numbers of A (such as momentum). The supertrace includes 
a relative minus sign for the fermionic component of the matrix, and we have explicitly 
displayed the trace over the discrete Matsubara bose [iv n = (2n)nT) and fermi (iu n = 
(2n + l)nT) frequencies. In the above expression, the self energy matrix S = Qq 1 — Q^ 1 . 
The quantity Y[Q] is the sum of all closed, two-particle irreducible skeleton Feynman graphs 
for the Free energy. 

As mentioned earlier, the Luttinger Ward functional Y[Q] has the property that its vari- 
ation with respect to Q generates the self energy matrix, 

SY\Q\ = -TStr[£5S]. (31) 
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The variation of the first term in F with respect to Q is given by 



5 (TStr [In (-^^ 1 )]) = TStr [-G^SG] 



(32) 



By using the relation £[£/] = g o — G , the variation in the second term in F is given by 



5 (TStr [([Go 1 - g-^g)]) = TStr [g^sg] 



(33) 



so that the total variation of F with respect to g 

s 



= TStr 



g- 1 + g Q l -e 1 56? 



identically vanishes, 



50 



0. 



(34) 



(35) 



Now the Hubbard Stratonovich transformation that we have carried out on the interaction 
Hi assures that the only place that the coupling constant olJk enters, is in = — -^j-.- 

This means that in the r.h.s of Eq. ()29j) a enters explicitly only through S x = — (a Jr) 1 — 
[g x ]-\ Then, 



dF 

da 



SF dg dF 
5g da da 
=o 

T 



(a 2 J, 



-Tr 



K 



g x {q,iuj n )e 1 ' 



(36) 



But by comparison with (|28j) . we see that = anc l since, F(a = 0) = F(a = 0) in the 
non-interacting case, the two quantities must be equal for all a, i.e 



F = TStr [In (-gT 1 ) + EG] + Y[g}. 



There are various points to make about this derivation: 



(37) 



The derivation is very general. Its correctness only depends on the stationarity of 
F\g] with respect to variations in g and the equivalence between dF/da and dF /da. 
This means that the above expression will hold for broken symmetry or "Higg's phase" 
solutions that involve off-diagonal components to £7[20|. In the context of the Kondo 
lattice, this means that g can be extended to include anomalous boson pairing terms 



13 



that are driven by short-range antiferromagnetic correlations, or alternatively, off- 
diagonal terms driven by long-range magnetism, which mix the conduction and \ 
fermions. This expression can also be used to describe superconducting states, where 
the conduction electron and % propagators will contain off-diagonal terms. 

Even though the free energy is a functional of three Green's functions, there is no 
overcounting. 

The above Free energy functional can be used as a basis for developing conserving ap- 
proximations that generalize the Kadanoff Baym approach to a constrained system[25]. 
In particular, in the large TV limit, the skeleton graph expansion for Y truncates at 
the leading diagram (Fig. providing the basis for a controlled treatment of both 



the magnetic and the paramagnetic phases of the Kondo lattice 



28). 



III. ENTROPY FORMULA 



In this section, we derive an approximate formula for the entropy of a Luttinger Ward 
system which becomes exact in the large N limit, and in any approximation where the vertex 
corrections can be neglected. The result, which we shall derive below, is 



S(T, Q) 



du. 



-Tr 



7T 



B,F 



dn(uj) 
dT 



(Imln [-£ _1 (u; - id)] + Im£(w)Re0(w)) 



where 



n{uo) 



n[uo\ 



(38) 



(39) 



is the matrix containing the Bose or Fermi-Dirac distribution functions n(u) = [e^ u — 
f(uj) = [e@ w + l] -1 . The trace is carried out over the Bose and Fermi components of both 
expressions. Here Q{oS) and £(&>) are the self-energies that have been analytically continued 
onto the real axis. This expression is extremely useful, for it only involves the low -energy 
part of the Green's functions of the particles, and does not involve the functional Y[Q]. A 
version of this formula was first quoted in the context of lattice gauge theory by Blaizot et 
al.Q. 

Key to our approach, is the notion that the Luttinger Ward functional F[Q] can be 
rewritten in terms of the real frequency Green's functions. This has the advantage that one 
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does not have to deal with the temperature dependence of the Matsubara frequencies. The 



approac 
systems 



r that we adopt is in fact, reminiscent of the Keldysh approach for non-equilibrium 



34 



35j . To preserve the symmetry between bosons and fermions, it proves useful to 



introduce a kind of "Keldysh" notation, writing 



1 



7]n{uj) 



(40) 



2 - m 

Use of this function preserves the "supertrace" symmetry of our approach. Consider the 
first part of the Free energy 



F 1 = TStr [In (-^T 1 ) + E0] . 



(41) 



We replace the summation in this term by an integral over n(u) around the imaginary axis, 

—ImStv[h(u)A(u - id)] (42) 

where Str [A] = TrL4e] — Tr[v4^] is the supertrace over the spatial and internal quantum 
numbers of A. By distorting the contour around the real axis, we obtain 

duj. 



■Im Str 



7T 



h(u>) {in {-g- 1 ) + {g Q l - g- l )g} t 



If we vary the Green's function, then we obtain 

duo 



5 Fx = - 



-ImStr 



7T 



(43) 



(44) 



where we have employed the notation S = ^ 1 — Q 1 In a similar fashion, when we vary g 
inside Y[Q], we obtain 



SY I — Im Str 

7T 



h(u) 6G(u) 



(45) 



The condition that the two terms cancel sets g^ 1 — g^ 1 = E, defining both the real and the 
imaginary parts of the self-energies in terms of the Green's functions. 

We can exploit the stationarity 5F/5g = to simplify the differentiation of the Free 
energy with respect to the temperature. When we differentiate F, we can neglect the tem- 
perature dependence of the spectral functions. 

dF OF 



S(T) = - 



dT 
15 



dT 



(46) 



The contribution to the entropy from F\ is then 

S x = J ^lmSti(^^[\n(-g-\uj-tS))+J:(u- l 5)g(u- l 5)]^. (47) 

The temperature derivative of the second term Y*[C/] in (|S7jl requires more careful consider- 
ation. The general diagrammatic contribution to Y contains a certain number of frequency 
summations. When we analytically continue these frequency summations, distorting the 
contour integrals around the branch-cuts of the Green's functions, we pick up "on-shell" 
contributions from each branch cut of the form 

h B {u)G" B {u) (48) 

and 



h F (u)G" F (u). (49) 

not hing 



These combinations are not hing more than the "Keldysh" Green's function well known from 
non-equilibrium physics (s^. \?>\ . In equilibrium, the Keldysh Green's function satisfies the 
fluctuation dissipation theorem 



G k {uj) = [G R (u) - G A {u))2h{u) = -AiG"{u)h{u) 

where Gr and Ga are the retarded and advanced propagators, respectively. In other words, 
hQ"{uS) = \Gk(u). The point is, that we can formally imagine evaluating Y in a general 
steady-state, replacing Gk — ► —Aig" {u)h{uj) at the end of the calculation. It is however, 
equally consistent to compute Y in equilibrium, and substitute g"{uj)h{uj) — > g K . The 
functional Y is closely related to the expectation value of the interaction energy. If we 
associate an amplitude \fa with the strength of each vertex, then the expectation value of 
the interaction energy has the expansion 

Qr, Ga) 

where Qk,r,a are the fully renormalized non-equilibrium propagators. The functional Y is 
the weighted sum 

y[q k , g R , g A ] = V — Y n {g K , g R , g A ) 

— ' n 

In this way, we can absorb all thermal functions into the Green's functions, so that all 
temperature dependence is entirely contained within the Keldysh Green's functions. To 
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differentiate Y with respect to temperature, we need to determine 5Y/5Qk- Unfortunately, 
Y is not a true generating functional for the Keldysh self-energies, and we must be very 
careful in taking the next step. 

Once we know the variation of Y with respect to the Keldysh Green's function, we can 
immediately compute the variation with respect to temperature. Now, from (J45|) . we see 
that we may write 



5Y 

If we now identify 



-Str 



7T 



J?(u)h(u)6g"(u) + h{Lo)Y!'{Lo)5Q'{Lo) 



(50) 



h(u)g" 



IK 



(51) 



as the Keldysh components of the propagator and self-energy, respectively, it is very tempting 
to write 



5Y = I ^Str 



7T 



(52) 



This form would always be correct if Y were a true generating functional for Keldysh propa- 
gators. Unfortunately, in the Keldysh approach, the distinction between the "measurement" 
and "response" vertices means that in general, the differential of Y with respect to Qk does 
not generate the real part of the self energy. Fortunately, this difficulty vanishes in the 
leading large N approximation, which is sufficient to generate a wide class of non-crossing 
approximation schemes. Diagrams involving vertex corrections do not satisfy this relation- 
ship, and there are corrections to the above diagram. Interestingly enough however, our 
large N expansion does absorb all the "RPA" diagrams of the interaction lines into explicit 
propagators, enabling the entropy formula to be derived for interacting collective modes and 
fermions. 

With this provise, we write 



5Y 



4 v ' 



where E' = \ - i6) + E(cj + iS)). 



It is worth making a short diversion at this point to explicitly demonstrate that this 
relationship works in the large N limit. Consider the leading order form for Y in the large 
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N expansion, 

2 



Yi = K x N x T 2 G c (k)G x (q - k)G b {q) (53) 

^ k,q 

where k = (iu n k), q = [iu n ,q) is a shortened space-time notation. When we carry out the 
Matsubara sums in this expression, we obtain 

Y = -K j°° ^ Yl \G' c G'^Gl{h x h b - 1) + G"G' x G b h c h b + G^G' 6 /a x } , (54) 

where we have used the short-hand notation G c = G c (k), h c = h c (u), G x = G x (g — fc) etc, 
which can be rewritten as 

^ = — / — — Re {G r c G k x G Kb + G KcG r x G Kb + G k c G k x G m + Gij x GR C GAb} , 

fc,<f 

(55) 

This last expression can also be derived using Keldysh formalism to calculate the leading 
order expression for (Ht) jsq. The variation of this expression with respect to Gk gives the 
leading order expressions for the boson, \ an d conduction self-energy, thus for instance, 

5Y i , f du) r 



[6G Kb {q) 4 



= E / - - fc ) + G^(*)fc(*)Gi( g - fc)] (56) 

a- 

where the replacement K — > /c = occurs because K contains the contribution of all iV spin 
channels of G^, and we are only differentiating with respect to one of them. This result 
can be independently confirmed using Matsubara techniques. 

Let us now return from this digression and continue to calculate the temperature depen- 
dence of Y . When we differentiate the general expression with respect to the temperature, 
we only need to keep track of how each of the thermal functions changes. This leads to the 
contribution, 

h{u)g"{u) = % -Q K {u) -> h{u)Q"{u)) + ^-G"(u)5T 

= % -[Q K (u)+8Q K (u)) (57) 



where 



{sGkH = d -^G"{u)8T (58) 
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We can now combine these results to obtain 



dY 
df 



duj 

7T 



Str 
Str 




(59) 



so that 



dY 



^Str 

7T 



dh(u) 
~~dT 



ReL{u)G"{uo) 



(60) 



This is a key element in our proof of the entropy equation. 

When we add S\ and S 2 together, S 2 partially cancels the second-term in Si, yielding 
the final answer 



S — Si + S2 



du , 



-Str 



7T 



dh{u) 
~~dT 



(Imln [-g-\z)] + ImE(z)Re£(z)). 



(61) 



We may replace h — > n and Str — > Tr b,f to recover (jHHjl . There are a few important points 
to be made about this result 

• Our ability to write the entropy in terms of the Green's functions only works because 
in the leading order approximation for Y, there are no hidden collective modes which 
carry the entropy. 

• This approach is general, and can, for example be applied to the interacting electron 
gas, and various interacting plasmas, such as the quark gluon plasma[21]. By treating 
the interaction line as an independent particle (photon) , the leading order approxima- 
tion generates a generalized RPA-Eliashberg scheme for the self-consistent interaction 
and electron self-energies. Note that the entropy formula can not be used in an RPA 
scheme in which the interaction line is not treated as an independent particle 22], and 
it does not work if one includes vertex corrections. 

• It may be possible in future work to evaluate a more general expression for SY/ 5Qk{^), 
permitting one to generalize the entropy formula to higher order approximations for 
Y. 
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IV. CONSERVATION LAWS AND SUM RULES 



Sum rules are intimately related to the existence of conserved charges. For example, in 
the Kondo lattice, there are three independent types of conserved charge: the electric charge 
(in each channel), the total spin and the conserved number of bosons at each site. Each 
conserved charge corresponds to a gauge invariance of the Lagrangian. For example, the 
electric charge (for channel v) is 

^ = E4a ( 62 ) 

k,a 

and this quantity is associated with the global gauge transformation 

H v « - e** (T) <w Xiu - e~^) X]v (63) 

The opposite sign in the exponents expresses the opposite charge of the two fields. Similarly, 
the conserved number of bosons at each site 

= ( 64 ) 

a 

generates the local gauge transformation 

In phases where these symmetries are unbroken, they give rise to sum rules, which we now 
derive. 

The general form of these gauge transformations is given by 

V> c -> e^ r ^Vc ( 66 ) 

We can relate this gauge invariance to the conserved charge by examining how the time- 
derivative terms in the action behave under this transformation. In general, the time deriva- 
tive part of the action has the form 

S = / dr^jdr^ (67) 







where the elements 7^ of the diagonal matrix 7 are unity, 7^ = 1 for conventional fermions or 
bosons and zero, 7^ = for Hubbard Stratonovich fields without any short-time dynamics, 
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such as the Grassmanian field Xin introduced in the Kondo lattice. We can relate 7 to the 
frequency dependence of Qq 1 , 

7 = dg °" lM - (68) 
— duj 

Under a time-dependent gauge transformation, the change in the action is given by 



AS = ^ jf dT^e- i0 ^jd T e m ^ - S 

= i / dTd T 6(T)Q(T) 

Jo 



= -1 I dr9(r)d T Q(r) (69) 
Jo 

where we have commuted the diagonal matrices q and 7 to obtain 

Q = ii\qj)i>. (70) 

Invariance of the action under time-dependent gauge transformations, implies that d T Q(r) = 
0, i.e. the charge is conserved. In this way, we see that 57 is the single particle operator 
associated with the charge Q. Notice that the full charge operator only depends on those 
fields with a time derivative. When the system is probed at short times, only these fields 
carry the charge. However, at long times, it is the gauge charge q which determines the 
charge of physical low energy excitations. 

Using (jZOJ), we can now write the expectation value of the the charge Q in the following 
compact form 

(Q) = (iiMa) = -^Str [qj_g] , (71) 

At zero temperature, provided there is a scale to the excitations, then the summations over 
discrete Matsubara frequencies can be replaced by a continuous integral. In particular, 

— Sti[A(cu)]e^ + (72) 

-too 2VT2 

where Str[A] without an underline indicates a supertrace purely over the spatial variables 
of A. This enables us to write 

/ioo j 
— Str [TfGiu)] , (73) 

where for clarity, we have temporarily suppressed the convergence factor e^ 0+ . Using flUD, 
7 = {dQo )/{du), we then obtain at T = 



<o> = - / ^ 
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' JGq\ 
dui 



(74) 



Now let us consider the effect of the gauge transformation on the Free energy. If we carry 



has the effect of shifting the frequencies associated with each propagator, oj^ — > + q^Auj, 
where is the gauge charge of the propagating particle. In the zero temperature limit, this 
shift in frequency can be made infinitesimally small, Au^ — > Sui^. 

Diagrammatically, the conservation of charges corresponds to the existence of closed loops 
(Fig. El) in the skeleton digrams of Y through which the conserved charge circulates. The 
gauge transformation causes the frequency variable in each of these closed loops to shift by 
an amount qAuo. This change does not affect the value of Y, and vanishing of the change in 
Y then gives rise to new Ward identities. 



FIG. 3: Illustrating the shift in frequencies around a closed loop in Y[Q] for the Kondo lattice, 
generated by the gauge transformation associated with charge conservation. 

When the frequency around the closed loops is shifted, each propagator entering into Y 
changes, with 8Q^{uj) = ^q^5u, so that the change in Y at zero temperature is 



out a gauge transformation ip£ — > e^^^tp, 





/ 





(75) 



so that 
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Integrating by parts, we obtain the Ward Identity 



3 du 

Combining eqs. (J73J) and ()76|) we obtain 



-Str 



du 



(76) 



(Q) 



-Str 



2vn 
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.0+ 



(Str [glB^" 1 )])^ 



,0+ 



(77) 



where we have explicitly reinstated the convergence factor. Folding the integration contour 
around the negative real axis, (as dictated by the convergence factor) we then obtain 
"° du d 1 



(Q) 



7r du 2i 



(Str [qlni-g-^u-id))] - Str [gh^-fiT 1 ^ + i8))]) 



/0 J, . J 
— -f-Im Str \q\n(-g- l (u-i8))\ 
-oo 7T du L 





7T 



Im Str [glnC-^-^w-itJ))] 



(78) 



-oo 



-Im Str [gln(-^ 1 (cj-^))] 



(79) 



w=0 



The presence of the minus sign inside the logarithms here is chosen so that the lower bound 
of the integral vanishes. We are then able to write 

7T 

As it stands, the trace in this expression involves both bosonic and fermionic components. 
The latter will give rise to Fermi surface volume contributions in the lattice, and to phase 
shift terms in the case of impurity models. What about the bosons? General arguments lead 
us to believe that the bosonic component in the trace vanishes. The important point here, 
is that the right-hand side of this sum counts the number of one-particle states that drop to 
negative energy. Whereas fermions can acquire a negative energy through the formation of 
a resonance or Fermi surface, bosons can not - they only condense. If the bosons condense, 
then the symmetry they are associated with is broken, and we can no longer apply the sum 
rule. If they do not condense, then we expect the associated phase shifts to be zero. The 
final form of the sum rule is then 



(Q) = -Im Tr 

7T 



q\n(-G F 



(80) 



=0— i<5 



where for clarity, we have suppressed the frequency argument of the Green's function. 
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V. KONDO IMPURITY MODEL 



As a first illustration of the sum rules, let us apply them to the Kondo impurity model. 
There are two sum rules to consider here corresponding to conservation of Q c and n b = 25 = 
K in the case of the perfectly screened solution. Let us first consider 



n b 



(81) 



Even though the x fermions do not enter into this conserved charge, they carry the gauge 
charge, and we have (q c , q b , q x ) = (0, 1, 1). Applying the sum rule (J7HJ), we deduce 



{n b ) 



-— Im Tr 

7T 



H-S-, 



H — Im Tr 



7T 



(82) 



where for the moment, we have retained the bosonic trace. When we carry out the trace 
over the spin components of the boson and the channel components of the x fermion, we 
obtain 

(n b ) = -N^ + K^, (83) 

7T 7T 

where 



5 b = Im [ln(A + E 6 (0 - iS) + i5)} 

<\ = Im ln(^ + S x (0-z<5)) 
Jk 



(84) 



are the "phase shifts" associated with the boson and holon fields. Note in passing that at 
zero temperature, E&(0 — iS) and S x (0 — iS) are purely real. The argument of the logarithm 
in 5b, is the renormalized chemical potential A* = A + £&(0) of the boson field. In the "filled 
shell" configuration 2S = K, the formation of a perfectly screened Kondo singlet generates 
a gap for the addition of extra Schwinger bosons to the ground-state, which forces A* > 0, 
and 5b = gued previously. It follows that 

.5, 



K—. 

TT 



(85) 



But n b = K is part of the constraint, so it follows that 5 X = n. This result can be understood 
by relating 5 X to the sign of the renormalized Kondo coupling constant, as follows 



— = -Im 

7T 7T 



lnC-L + E^O-i*)) 

"K 



9(-J* K ). 



(86) 
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In this way, we see that 5 X = tx corresponds to a residual "ferromagnetic" coupling. This 
is consistent with our expectations, for when an additional Schwinger boson is added to 
the Fermi liquid, it increases the impurity spin by one half unit to form an underscreened 
Kondo model, where the residual spin coupling is indeed ferromagnetic. By contrast, in an 
antiferromagnet or spin liquid, the local moments spins pair-condense, mutually screening 
one-another. In this situation, an additional Schwinger boson at any one-site, will now be 
free to undergo a Kondo effect with the conduction electrons, so the residual spin-coupling 
is antiferromagnetic. These results are also confirmed in the large N limit. 

Let us now turn to the sum rule derived from conservation of electron charge. Here the 
conserved quantity is 

fl.=EC^ (87) 

and the corresponding gauge charges are (q c , qb, q x ) = (1, 0, —1), so the sum rule becomes 



Q e = i (im Tr [ln(-^ 1 )] - hm Tr [ln(-^ 1 )]J 



(88) 



=0—i<5 



Now in an impurity model, it is convenient to subtract off the total charge in the absence 
of the impurity, given by 



Q( 0) = -Im Tr 

7T 



H-e, 



cO / 



(89) 



=0-i<5 



where Q c0 = is the non-interacting propagator of the conduction electrons. If we subtract 



this from (|88j). we can combine the conduction electron traces and replace ln(— Q c x ) — 
ln(— GcQ 1 ) = ^(GcoQc 1 ) = l n (l — Gco^c)- The change in total charge is then given by 



Qe-Q 



(0) 



-Im Tr 



TX 



ln(l - GcqE c 



1 



Im Tr 



TX 



(90) 



According to the Anderson-Clogston compensation theorem, the total change in electron 
charge due to a Kondo or Anderson impurity vanishes in the infinite band- width limit, 
and at finite band-width the change is of order the ratio of the Kondo temperature to the 
bandwidth AQ C = 0(Tk/D), and can be neglectedQ, Q- Setting AQ e = Q e - Q° e = 0, 



we obtain 



where we have identified 



AQ e = 



5 C = Im Trg 



5 r 5 V 

NK- - K— 



ln(l - g c0 T, c 



(91) 



(92) 



-id 
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where the trace is purely over momenta. If we expand the logarithm using a power se- 
ries, we see that the trace over internal momenta in (}9*2*)) is accomplished by replacing the 
momentum dependent Green's function Gc (k, u) by the local electron Green's function 

0o M = £fc^(fc> w )> so that 



S c = Imln [1 — gQ(u)H(u)] 



(93) 

=0-i5 



In the large band- width limit, go(u — iS) — > inp, where p is the density of states per spin 
per channel, so that 

5 C = Imln[l -mpE(0)) (94) 

There is thus a direct link between the phase shift of the \ fermion and that of the 
conduction electrons. Combining this with the earlier result 8 X = n, we obtain 

5 C = ^ = * (95) 

Notice that 

• in the large N limit the phase shift identity permits us to relate the conduction electron 
phase shift 5 C , which is 0(1/N) to the \ phase shift 8 X , which is finite in the large N 
limit. In this way, the sum rules enable us to study the phase shift and Fermi surface 
volume changes of the Kondo impurity and lattice, in the large N limit. 

• the development of the conduction electron phase shift does not occur because conduc- 
tion electron states drop beneath the Fermi sea: it is a consequence of the injection of 
new quasiparticle states into the Fermi sea as a response to the formation of spinless, 
charged holons. 



VI. SUM RULE FOR THE KONDO LATTICE 



When we come to consider the Kondo lattice, there are a number of additional subtleties 
that must be entertained. One of the most important aspects of the discussion here concerns 
the possibility that the extremal solutions to F[Q] break the U(l) symmetry of the bosons, 
to produce "Higg's" phase in which the Schwinger bosons either pair condense, or hybridize 
between sites. In the large N solutions to the model, these phases appear to develop as a 
precursor to the formation of the antiferromagnet, where the bosons themselves condense. 
The main effect on the sum rules is two fold: 
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• we can not assume that the x fermion is localized. This implies that we must consider 
the momentum dependence of the x fermion. 

• once the Schwinger bosons pair condense, the Ward identities associated with the 
conservation of rib no longer apply, because the Luttinger Ward functional Y[Q] is no 
longer invariant under shifts of the boson frequency. 

However, we can still take advantage of the Ward identity associated with charge conser- 
vation. Provided that the Schwinger boson is uncondensed, i.e, there is no magnetism, then 
charge conservation guarantees that 



duo 
2tH 



Tr, 



—Qq 



0, 



(q c = l,q x = -1) 



(96) 



Now in the paramagnet, the ho Ion and conduction electron fields are unmixed, so we can 
replace 

Tr F -> Tr c + Tr x . 

This is a special property of the paramagnet. By contrast, in the antiferro magnet, the 
Schwinger bosons are condensed, so the holon and conduction electron fields admix, and 
both X and Q contain off-diagonal terms so that Trp remains a single, integral trace over 
the two admixed fields. 

For the paramagnet, we can write 



= 



dcu 



2ni 



Str 



dcu 



Qq 



duo 
2^1 



-Tr 



dui 



du 



2ni 



Tr 



r°° duj_ r d D k 

J_ ioo 2riJ (2^ 



N 



Xn 

dE c (k, oj) n N dE x (k,uj) 



duj 



-G c {k,uj) - 



duj 



G x {k,u) 



(97) 



It is this identity that permits us to extend the Luttinger sum rule to the case of the 
paramagnetic Kondo lattice. In the paramagnetic Kondo lattice, the charge sum rule 



Qe = ~ (Im Tr [ln(-^ 1 )] - Im Tr [ln(-^ 1 )] ) | w=0 _ 



(98) 



now involves a trace over momentum: 



% = N V -Im In [e n + E c (jfe, - i6) + iS] - V -Im In + E x (k, - iS)] 



(99) 
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The first term in this expression is the electron Fermi surface volume, 

E ^Imln[e fe - + X c (k,0 - i6) + i5] = £ B{-E % ) = (100) 
k k 

— * 

where the region where E% = + ReS c (/c, E£j is negative defines the interior of the Fermi 
surface. The second term can be interpreted in a similar way - the momentum trace over 
the logarithm of Q' 1 

E ^ Im ln [ + s x(^ - tf)] = E d i- J K(k)} = j^jp (101) 

k k 

can be seen as the volume of the region in momentum space where the effective interaction 
J^{k) = —Q x (k,0)~ 1 is negative, or ferromagnetic. Following our earlier discussion, the 
ferromagnetic sign of the residual interaction is a consequence of the fact that additional 
spins added to the state completely decouple from the Fermi sea. In the simplest scenario, 
J^(k) < for all k, in which case — L an d the sum rule becomes 

-t^- 1 < 102 > 

where n e is the electron density per unit cell, per conduction electron channel. It follows 
that the total Fermi surface volume expands by one unit per unit cell, 

= n e + l (103) 

Now if there is some region of momentum space where J^{k) is not ferromagnetic, it 
follows that the x fermions will have a Fermi surface, and the excitation spectrum will now 
involve charged, spinless fermions. In this phase, 



K (2ir) D (2vr) J 



= - = " 7^75 ( 104 ) 



Such a "spin-charge decoupled" phase could not develop from a Fermi liquid without a phase 
transition. So long as the heavy electron fluid does not undergo such a phase transition, the 
Fermi surface must remain large. 

Let us now consider what happens in the antiferromagnet. In the general higher dimen- 
sional case, the boson field condenses, causing the conduction and x fields to hybridize to 
produce a single species. In the antiferromagnet, it is more logical to make a particle- hole 
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transformation of the x field, writing ip^ = x\ £■ 111 this case, = q c have the same charge, 
and the sum rule becomes simply 



Q e = hmTr F [ln(-^ 1 )] 



(105) 



where Qp is the admixed propagator for the combined conduction and ip fields. The right- 
hand side can not be separated into x an d conduction parts, and as such, defines an admixed 
set of Fermi surfaces, with an average Fermi surface volume which counts the total charge 
per unit cell 

n, = N^0. (106) 

The Fermi surface volume for the antiferromagnet and the paramagnet must then differ by 
1/N th of a unit cell. In the special case of N = 2, if the spatial unit cell doubles, this formal 
change in Fermi surface volume is indistinguishable from the new size of the Brillouin zone. 
However, even in this case, we can imagine more general classes of antiferromagnet, such as 
an incommensurate helimagnet, where the unit cell size is unchanged, but the total Fermi 
surface volume must jump. 



VII. SPIN-CHARGE SEPARATION AND QUANTUM CRITIC ALITY. 

Although the sum rules do not provide us with any details of the dynamics, they provide 
stringent constraints on the way the spectrum of excitations can evolve in the Kondo lattice 
as we approach the magnetic quantum critical point. 

One of the fascinating aspects of our sum rule 

= cm* (107) 

is that it suggests that a Kondo lattice may develop low-lying, spinless charged fermionic 
excitations near a quantum critical point. This is clearly a controversial idea, the conse- 
quences of which we now explore in this discussion. If we take this idea seriously, then two 
important points seem to emerge: 

• In the heavy fermi liquid, spinons have "ionized" into electrons and a background filled 
sea of holons. 

b a ^e-+ X + (108) 
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It is this ionization process that lies behind the expansion of the Fermi surface and 
the sum rule. 



• There is an intimate link between the Kondo interaction and the formation of holons. 
The link between the propagator of the holons and the renormalized Kondo interaction 

-QM w) = i * n v = J*k(Q, w) (109) 

implies that poles in the holon spectrum correspond to divergences in the Kondo 
interaction. As magnetism drives the Kondo interaction back to weak coupling, we 
expect J*(q,u) to pass through infinity, so the holons should become gapless at some 
point in this process. 

We can in fact use the Luttinger Ward approach to gain insight into the low-energy 
thermodynamics of the Kondo lattice Q| . Provided the replacement of frequency sums by 
continuous integrals is valid at low temperatures, Luttinger argued that the low energy 
thermodyanmics is determined by the leading logarithm in the Free energy functional 
evaluated with the zero-temperature Green's function. In our case, this leads to the relation 

F(T)^^TStr[ln(-^ ;)], (T - 0) 

where is the zero temperature Green's function. Carrying out the frequency integrals, 
this implies that the low temperature entropy is given by 

spinon and holon entropy 

Fermi liquid 



r d D k r du ( df „ 1/r A 

where the zero temperature propagators are to be used. The last term in this expression 
provides the T— linear entropy of the Fermi liquid, the first two terms are the "spinon" and 
"holon" contributions. For a Landau Fermi liquid to form, these last terms must cle arly 
be gapped at low temperatures. This feature is observed in the large N solutions 28] . 
What it also makes clear however, is that the holons can not become gapless without a 
thermodynamic departure from Fermi liquid behavior, which implies some kind of quantum 
phase transition. 
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However, even though the holons are gapped in the Fermi liquid, they can exist as low- 
lying spinless, charged excitations. It is particularly interesting to speculate that the gap to 
holon formation closes at a quantum critical point. This is a possibility that can be studied 
in the large N limit, but the general arguments should hold at finite N, so long as the 
Luttinger Ward approach is valid. One way to go beyond the large N limit, is to to carry 
out a fully self-consistent treatment of the leading Luttinger Ward Free energy functional, 
updating the conduction electron self-energy and feeding the full conduction electron Green's 
functions into the calculation of the self energy for the x an d Schwinger boson field. Such 
an approximation goes selectively beyond the large N limit, but it will satisfy the Ward 
identities. Moreover, assuming that the "filled shell" stability of the fully screened ground- 
state is not an artifact of the large N limit, we also expect the gap in the Schwinger boson 
and the \ fermion particle spectrum to be preserved at finite N. 

Another way to check the picture emerging from this approach would be to examinelow- 
dimensional systems. Various authors have examined the possibility of a transition from 
small to large Fermi surface in the one-dimensional spin- 1/2 Kondo latticeQ, Q- The 
difficulty with this model, is that the gaplessness of the Heisenberg spin-1/2 chain tends to 
make the Kondo coupling relevant, no matter how small the Kondo temperature. A model 
that avoids these difficulties is the 5 = 1, two channel Kondo lattice. This model involves a 
Haldane spin-1 chain coupled via Kondo interactions to two one dimensional, non-interacting 
Hubbard models, as shown in Fig. 0] 

When the Kondo interaction is weak, the Haldane chain is gapped, and the "Fermi 
surface" of the Kondo lattice is small, giving rise to Friedel oscillations at wavevector Q = 
2k p\ When the Kondo interaction is strong, a Kondo lattice with a large Fermi surface is 
expected to form. If the "large" and "small" Fermi surface phases, where Q = 2kp + n. 
The way these two phases are linked is particularly interesting. If they are linked by a single 
quantum critical point, then we expect Q to simply jump at this point (Fig. HJd.). However, 
if the holons form a gapless phase then over an intermediate range of parameters, we would 
expect the Friedel oscillations to exhibit an intermediate wavevector (Fig. Et-), 

Q = 2k° F + xn. (Ill) 
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v=2 



(b) V FS 2kP 2kP+n 

q T K / J H 

(c) V FS 2kf ^P+xn 2k p + iT 

I W////////M* 

T K /J H 

FIG. 4: (a) The two-channel, spin 5=1 Kondo chain, (b) In the simplest scenario, the transition 
from a Haldane chain with small Fermi surface, to the heavy electron phase with large Fermi 
surface will occur via a single phase transition, (c) The sum rules allow for an intermediate phase 
with a Fermi surface that lies between these two extremes. 
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